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Abstract. Recently two anisotropic superconducting gap functions have been observed in the skutterudite 
PrOs4Sbi2. These order parameters are spin-triplet. There are at least 2 distinct phases in a magnetic field, 
bearing some resemblance to superfluid 3 He. Here we present an analysis of the thermodynamic properties 
in these two superconducting states within the weak-coupling BCS theory. 

PACS. 74.20.Fg BCS theory and its development 



1 Introduction 

Superconductivity in the body-centered cubic heavy-fermion 
(HF) skutterudite PrOs4Sbi2 was discovered in 2002 by 
Bauer et al |H2l3| . Since then many experimental and the- 
oretical studies of this compound have been reported. This 
compound possesses several interesting and unusual char- 
acteristics: two distinct phases (the A phase and B phase) 
in a magnetic field, nodal superconductivity with point 
nodes, and triplet pairing with chiral symmetry breaking 
|6I7I8| . The phase diagram is still controversial. In Fig. 1 
recent measurements by Measson et al [2j are shown. 

It was recently observed that the magnetothermal con- 
ductivity data | 4I7| in this compound are consistent with 
anisotropic superconductivity using the gap functions 
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Fig. 1. Phase diagram by Measson et al|3j. The top line and 
points are the upper critical field H C 2, while the lower ones are 
the phase boundary H'. 



(k z +ik x )/Jkl + kle^ = {k x +ik y )/Jkl + kl. 



The 



factor of 3/2 ensures proper normalization of the angular 
dependence of the order parameter. In Eq.(2) the nodal 
direction is chosen to be parallel to (001). 

We note that the proposed order parameter (2) lies 
outside of the usual classification schemef!!:, in which or- 
der parameters correspond to a single irreducible represen- 
tation of the rotation group. However, this hybrid order 
parameter appears to be necessary to reproduce the ob- 
served B-phase gap structure. A similar situation has been 
observed in the borocarbide superconductors |12| . 

The cubic symmetry of PrOs4Sbi2 suggests order pa- 
rameters which are invariant under the T^ cubic tetra- 
hedral symmetry group applicable to this crystal as 
well as reflections (containing the origin) about the planes 
of the crystal parallel to the cube faces |T5|. As suggested 
in JOjj one possible invariant is k x k 2 + kyk z + k 2 z k x . This 
belongs to the Ai representation of T/j [IT)] . 

This combination can be recast as 1 — — ky — k\ , 
thus forming the basis of the proposed order parameter of 
the A phase. Furthermore, in weak-coupling BCS theory 
the quasiparticle density of states and the thermodynam- 
ics depend only on | /" | |X (Sj - the magnitude of the angle- 
dependent part of the order parameter. For this reason, 
an order parameter which breaks chiral symmetry still re- 
tains the essential features of the cubic symmetry, and 
is in fact necessitated by the triplet pairing observed in 
this compound. |fil7l8| Triplet pairing requires that the or- 
bital wavefunction be antisymmetric under particle inter- 
change. The phase factor proposed in Eq. (1) meets this 
requirement. 

The proposed B-phase order parameter breaks the cu- 
bic symmetry more manifestly. Nevertheless, there is al- 



2 



David Parker et al.: Anisotropic superconductivity in PrOs4Sbi2 





Fig. 2. A-phase (left) and B-phase (right) order parameters 



most certainly a relationship between the A and B phase, 
particularly since the zero-field transition temperatures 
are so nearly equal. The simplest relationship consistent 
with the nodes at [001] and [00-1], despite requiring a hy- 
brid representation, would suggest |/| ~ 1 — h*. This order 
parameter again would have a phase factor included in / to 
ensure antisymmetry under interchange. These proposed 
order parameters are illustrated in Fig. 2 |17| . 

We note that while the B-phase is the prevalent phase 
in zero magnetic field, the A-phase exists at all tempera- 
tures below T c for fields between H* (the phase boundary) 
and H c2 . 

Below is a comparison of the predicted B-phase DOS 
with scanning tunneling microscopy (STM) data taken by 
Suderow et al JTBJ at T=0.19 K. The predicted B-phase 
DOS differs somewhat from that presented in [7] due to 
the use of the angular-dependent quasiparticle density-of- 
states, as well as an accounting for the energy and direc- 
tional resolution of the STM . Here we have assumed that 
the STM performed measurements along the nodal direc- 
tions, where a is the size of the momentum cone within the 
STM's spatial resolution. We note that the nodal structure 
can easily be masked by performing STM along a limited 
number of directions of single crystals. Here we have as- 
sumed A to take the B-phase weak-coupling value of 3.3 
K. We note that the observed small DOS for E < A/3 
can be reproduced by choosing a to be 3.0. We observe 
fair agreement, with some differences apparent surround- 
ing the quasiparticle peak at E = A. In the following 
we analyze both phases over the entire temperature range 
from T = to T r . 




Fig. 3. Comparison of quasiparticle DOS in B phase with STM 
data of Suderow et al. |18| 



Then, within the weak-coupling theory the gap equa- 
tion is given by 



A" 1 = 2ttT < f >- x 



(3) 



fE r2 TP 



2 Weak-coupling BCS Theory 

Wc focus on the superconductivity in the A and B-phases 
of PrOs4Sbi2, using the A(k) given by Eqs. 1 and 2 with 
|d| = A(T). 



where A is a dimensionless coupling constant, Eo is the 
cut-off energy, and w„ is the Matsubara frequency. Here, 
for the A-phase / = |(1 — fc 4 — ky — A; 4 ), and (. . .) denotes 
J dfi/A-K. For the B-phase, / = 1 — z 4 , and (. . .) denotes 
f Q dz . . .. The frequency sum in Eq. 3 is cut off at cj„ = Eq. 
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Fig. 4. A(T) for the A and B phases 



In the vicinity of T = T c and T = OK, Eq. 3 and Eq. 
4 can be solved analytically. For T — > T c we obtain 



T c = — E e- 1/X = 1.136-B e _1/A 



A 2 {T) 



2</ 2 >(2^T c ) 2 ln(T/T c ) 



7C(3) < / 4 > 
where 7 = 1.78 .. . is the Euler constant and 

A{0)/T c = -exp[- < f 2 > x < / 2 ln(/) >] 

7 

= 2.364, A - phase 
= 1.938, B-phase 
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In the low-temperature regime T/A(0) <C 1 one obtains 

T 



\n(A(T)/A(0)) = - 7 -^( 
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In Fig. 4 numerical solutions of A(T)/A(0) are shown for 
both phases over the entire temperature range. 

The values of A(0)a and A(0)b obtained from the 
weak-coupling theory offer a possible explanation for the 
multiphase diagram of PrOs4Sbi2. The value of the con- 
densation energy at T = is given by 



E 



-- < |/| 2 > N A 2 



(12) 



where / is the angular-dependent part of the order pa- 
rameter and No the normal state density of states at the 
Fermi level. For the A-phase, < |/| 2 > = 3/7, whereas for 
the B-phase, < |/| 2 > = 32/45. This yields two distinct 
condensation energies: 



-i.rw (T c A ) 2 , 

T.32A^ (Tf) 2 . 



(13) 
(14) 



If one uses the experimental values [Hj for and as 
1.887 K and 1.716 K, respectively, one finds that E^ is 
slightly lower than Ejf , contrary to observation. On the 
other hand, if we assume T C A = T C B , which the simplest 
interaction would give, we find Eg < E$. The difference 
between our assumed and the measured T c would in this 
case be due to some unknown external perturbation, such 
as the effect of the crystalline electric field (CEF), not 
accounted for in this treatment. One possibility is that 
the CEF affects superconductivity in the two phases dif- 
ferently, resulting in a difference in measured transition 
temperatures. Indeed, if we assume that, absent such an 
effect, we would have = Tf we find Qg < f}£ for all 
temperatures. 

Upon evaluating A(T), the thermodynamics of the sys- 
tem can be analyzed following Ref.^Hl- Let us start with 
the entropy: 



S s 



dEN(E)(f]nf + (1 - /) ln(l - /)),(15) 



where / is the Fermi-Dirac distribution (1 + e 



0E\ 



with 



(3 = 1/ksT. N(E) is the quasiparticle density of states, 



N(E) = N Re( 



\E\ 



VE 2 - A*f- 



(16) 



The electronic specific heat can be derived from the en- 
tropy via 



C — T— 
Cs ~ 1 dT 



(17) 



In Fig. 5 we show C s /~fsT versus T/T c for both phases. 
Here 75 = 2ir 2 No/3 is the Sommerfeld constant. We find 
the jump AC fC at T c to be approximately 0.93 for the 
A-phase and 1.20 for the B-phase. Data from Vollmer et 
al[2I shows these jumps to both be of order one, so that 
our model is consistent with this data. In addition, as ex- 
pected, the low-temperature specific heat is predicted to 
be proportional to T 2 |S] for both phases. Unfortunately, 
the presence of a Schottky specific heat peak|5] at low 
temperature makes assessment of the T 2 prediction diffi- 
cult. 

Also the thermodynamical critical field H C (T) can be 
obtained from 



F S (T)-F N {T) 



dT(S s (T)-S N (T)) (18) 



"8^ 



(19) 



Here Ss(T) and Sat = jnT are the entropies in the 
superconducting and normal state respectively. We show 
D{Zr) = H c {T)/H c (0) - (1 - {T/T c ) 2 ) for both phases in 

Fig. 6. The function D(jr) for both the A-phase and B- 
phase cases is slightly larger than for the isotropic s-wave 
case and somewhat smaller than for the d-wave case|20|. 
with the B-phase case slightly larger than the A-phase 
case. 
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Fig. 5. Specific heat C a /-y s T for PrOs 4 Sbi 
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Fig. 6. Deviation of critical field from parabolic dependence 
for p+h, s, and d-wave superconductors 



While there are a few reports of H C 2(T) for the A- 
phase and H*(T), the phase boundary between the A 
phase and the B phase [rT4l9l21l22l23| . no experimental 
data are available for H C (T). 

Finally the superfluid density is given by 

Ps\\{0) Jo 2T 

E \ 



Re( z z 



y/E*-A*(T)p 



(20) 



and 



P*x(T) 
P. 



±(0) 2J 2T 1 1 ' 

Re((l~z 2 )- 



E 



y/E* - A*(T)f* 



(21) 



where Re(. . .) refers to the real part, and the subscripts 
|| and _L indicate parallel and perpendicular directions to 
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Fig. 7. Predicted superfluid densities for PrOs4Sbi2 
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Fig. 8. Superfluid densities for PrOs4Sbi2 



the nodal points. The superfluid density, as expected, is 
isotropic for the cubic symmetry-retaining A-phase, but 
rather anisotropic for the B-phase. These superfluid den- 
sities are shown in Fig. 7. In the low-temperature regime 
(T <C A) both Eq.(20) and Eq.(21) can be expanded as 



p s (0) 2 V A 

Pb\\(0) 

PsB±{T) 



4 '""'A 



P.x(0) 

Close to the transition temperature, we find 

Psa(0) 7</ 4 > V 1 cl 
= 1.393(-lnT/T c ) 

Ps\\(0) 21 -IV 1 c> 
= 0.9567(- In T/T c ) 



(22) 
(23) 
(24) 

(25) 
(26) 
(27) 
(28) 
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In the figure above we have compared — 



>nCO 



,(o) 



(29) 
(30) 

with the 

data taken from Chia et al[22], assuming that the nodal 
points in A(k) are aligned parallel to H. Rather satisfac- 
tory agreement is observed for T < T c /3. But the theoret- 
ical p s ±(T) vanishes linearly with T c — T in the vicinity 
of T = T Cl whereas Chia et al 24 found a p s ±(T) which 
vanishes with essentially infinite slope at T c . 

Recently Chia et al jJS] also reported magnetic pene- 
tration depth measurements for a range of dopings x from 
0.1 to 0.8 in the compound P^Osi-^Ru^Sb^. Over the 
range from x = 0.4 to x = 0.8, exponential temperature 
dependence of the superfluid density was found, indicating 
an isotropic s-wave gap function in this regime. Of direct 
interest for this work, the superfluid density was found to 
go to zero linearly for all dopings, with no hint of the es- 
sentially infinite slope found |24| in the pure case (x=0). 
In addition, the slope of these linear curves at T c does 
not increase dramatically from x=0.8 to x=0.1. Further 
experiments at doping ranges between x = and x = 0.1 
are highly desirable, to examine more closely the apparent 
transition from nodal to conventional superconductivity 
taking place in this system. It would also be of value to 
confirm the rather unusual "infinite-slope" behavior ob- 
served in the pure sample near T c . 



3 Concluding Remarks 

We have worked out the weak-coupling theory of the A 
and B phases of the heavy-fermion superconductor 
PrOs4Sbi2. A simple thermodynamic analysis offers an 
explanation for the appearance of the lower-symmetric B 
phase at lower temperatures. The present model leads to a 
fair description of STM data taken by Suderow et al |18| . 
In addition, the present model for the B-phase describes 
the superfluid density determined by Chia et al 24 for 
the low-temperature regime, if we assume that the nodal 
points in the B-phase follow the magnetic field direction 
in the field cooled situation 7 . Since the magnetic field 
is the only symmetry-breaking agent, this appears to be 
plausible. We will present the results of an analysis in the 
case of impurities shortly. 

We thank H. Won, P. Thalmeier, K. Izawa, Y. Matsuda 
and H. Tou for many useful discussions. S.H. acknowledges 
financial support through PRF grant 41757-AC10. 
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